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Potential non-relativistic QCD provides a convenient framework to study the three heavy-quark 
system. We have studied its symmetries under exchange of the heavy quarks and determined the 
leading ultrasoft contributions to the singlet static potential and to the singlet static energy. 
Specializing to an equilateral geometry, we have solved the renormalization group equations in 
order to resum the leading ultrasoft logarithms for the singlet static potential, which is thus known 
at NNLL accuracy. 
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1. Introduction 

We have studied the symmetries of the three heavy-quark system under exchange of the heavy- 
quark fields and their implications for the various matching coefficients, i.e. the potentials, of 
potential non-relativistic QCD (pNRQCD) for the three heavy-quark system. Moreover, we have 
calculated the ultrasoft corrections of order a s 4 In a s to the singlet static energy and of order In /I 
to the singlet static potential of a three heavy-quark bound state. Whereas this has been achieved 
for the case of QQ systems more than ten years ago [fin, the result for three heavy-quark systems 
are new [0J. 



2. Three heavy-quark composite fields 

Quarks transform under the fundamental representation, 3, of the (color) gauge group SU(3) C . 
According to 

3®3©3= 1©8©8©10, (2.1) 

the three heavy-quark product state, Qi(xi,t)Qj(x2,t)Qk(x3,t) (i,j,k= 1,2,3 denote color indices), 
can thus be decomposed into a singlet, two different octets and a decuplet with respect to SU(3) C 
gauge transformations at a common point R. Employing equal-time straight Wilson strings, 

(j)(y,x,t) = &>explig^ ds(y-x)-A(x+(y-x)j,o|, (2.2) 

where A = A a X a /2 is the color gauge field, X a are the Gell-Mann matrices, and 8? denotes path 
ordering of the color matrices, we write 

Qi (xi , t) Qj (x 2 , t) Qk (x 3 ,t) = fa (x u R,t)(t>jf(x 2 ,R,t) <fe< (x 3 , R, t ) 

( 8 

x^(x 1 ,x 2 ,x 3 ,0S,, /it , + £ Afl (xi,x 2> 13,00^4, 

I a=l 

8 10 >| 

+ £ Sa (x 1 ,x 2 ,x 3 ,?)0^ + £ A S (x u x 2 ,x 3 ,t)Af, fk , , (2.3) 
a=l 5=1 J 

where S ;7 „ 0^,0^ and ^fjk are orthogonal and normalized color tensors that satisfy the relations 

Sc 1 nAa*(\Ab sab r\Sa*r\Sb sab \S \S' sSS' 

-ijk&jk — 1 ' ±*ijk ^ijk — , KJ ijk KJ ijk — , ^jk^jk — , 

SijkQjfk = Sijk^jk = Sijk^fjk = QjfkQjjk = Qdjk^ijk = Qjj'k^jk = , (2.4) 

with a, b G {1,...,8}, and 8,8' G {1, , 10} [pp. For simplicity, we have omitted an explicit 

reference to R in the argument of the singlet, S, two octets, O 4 and s , and decuplet, A, fields. 
Besides the time coordinate t, we rather list the position coordinates (xi,x 2 ,x 3 ) of the heavy-quark 
fields in the order (from left to right) of their appearance on the right-hand side of Eq. (23). The 
same convention is used for the color indices (i,j,k). 

The heavy fields fulfill equal-time anticommutation relations, 

{Qi(x,t),Qj(y,t)}=0. (2.5) 
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Thus, it is natural to ask for the implications of Eq. (2.5) for the singlet, octet and decuplet fields 



in Eq. ( |2.3| ). In order to do this, we stick to a specific (matrix) representation of the color tensors, 
namely [^] 



iijk ~ ~^ e 'j k ' — ijk ~ 2 £ijl ^ 1 ' —<'j' k ~ ^/f ( £ i kl ^il + e 'kl^jl) ) 



(2.6) 



and 



Al _ A 4 _a10 
±±111 ~ ±± 2 22 — ±±333 



A 6 - 1 

- {123} ~ vi' 



,8 



±±{112} —±±{122} —±±{113} —±±{223} "±±{133} —±±{233} 



1 



(2.7) 



where the symbol {ijk} denotes all permutations of the indices ijk; all components not listed 
explicitly in Eq. ( |2.7[ ) are zero. Moreover, note that as all six possible orderings of the heavy-quark 
fields can be obtained by applying the following two independent operations on the product state 
Qi{x h t)Qj(x 2 ,t)Qk(x 3 ,t), 

(/): Q i (x l ,t)^Q j (x 2 ,t), (ii): Qi(x h t) -H- Q k {x 3 ,t) , (2.8) 

it suffices to determine the transformation behavior of S, Q^", Sa and A under (/) and (ii). We find 



(0 



5(xi,x 2 ,x 3 ,0 
A 5 (xi,x 2 ,x 3 ,?) 
Aa (xi,x 2 ,x 3 ,?) 
O 5fl (xi,x 2 ,x 3 ,0 



S(x 2 ,xi,x 3 ,?) 

-A 5 (X 2 ,X!,X 3 ,?) 

Aa (x 2 ,x 1 ,x 3 ,t) 
-0 Sa (x 2 ,x u x 3 ,t) 



(2.9) 



and 



(ii) 



S(xi,x 2 ,x 3 ,t) = 5(x 3 ,x 2 ,xi,0 
A 5 (xi,x 2 ,x 3 ,f) = -A 5 (x 3 ,x 2 ,xi,f) 

O Aa (x 1 ,x 2 ,x 3 ,0 = -cos(f)O Aa (x 3 ,x 2 ,xi,0 + sin(f)O Su (x 3 ,x 2 ,xi,0 
O & '( Xl ,x 2 ,x 3 ,0 = sin(f)O Afl (x 3 ,x 2 ,x 1 ,0 + cos(f)O ft '(x 3 ,x 2 ,x 1 ,0 



(2.10) 



3. The pNRQCD Lagrangian for the three heavy-quark system 

As pNRQCD for the three heavy-quark system is usually formulated in terms of singlet, octet 
and decuplet fields, the symmetry relations (19) and ( 2.10| ) have some immediate consequences 
for the form of the pNRQCD Lagrangian, which is organized as a double expansion in l/m and in 
the relative coordinates r, (i = 1,2,3). To zeroth order in both expansions, it reads [^| 

^PNRQCD = / d V <* 3 A { S ' [«b " V s ] S + A ' [ iD ° ~V d ]A + [/Do - Vf] & 
+O st [/Do - V|] s + c/ 1 " [-v&] o s + O st [-V£ s ] } 

+E«W ~jF^ V ■ (3-D 
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Equation (3.1) describes at zeroth order in the multipole expansion the propagation of light quarks 
and ultrasoft gluons as well as the temporal evolution of static quarks. Higher order terms, not 
displayed explicitly here, account for corrections due to finite heavy-quark masses and interactions 
between heavy-quarks and ultrasoft gluons. As there is no preferred ordering of the heavy-quarks, 



Eq. ( p.l[ ) has to be invariant under different orderings of the heavy-quark fields. 
The potentials V can be expressed in terms of the relative vectors 



xi - x 2 , 



r 2 = xi - x 3 . 



r 3 = x 2 - x 3 . 



(3.2) 



For the ordering of the heavy-quarks as in Eq. (^3|), i.e. S = S(xi,X2,X3,/), O 4 = A (x\,X2,x^,t), 
qS = s (x\ ,x 2 ,x 3 ,f) and A = A(xi,X2,X3,£), they are defined as V = V(r\ ,r 2 ,r 3 ). Given the sym- 
metry relations for the singlet, octet and decuplet fields, Eqs. ( |2.9| ) and ( [2.10 ), it is straightforward 



to also derive corresponding symmetry relations for the potentials in Eq. (3.1). The singlet and 
decuplet potentials remain invariant under (i) and (ii), whereas the octet potentials transform as 



V A °(ri,r 2 ,r 3 ) 
V s °(r h r 2 ,r 3 ) 
[VX s (ri,r 2 ,r 3 ) 



^(-r l5 r 3 ,r 2 ) 
^(- ri ,r 3 ,r 2 ) 
-V A ° s (-ri,r 3 ,r 2 



(3.3) 



and 



(ii) : 



/V£(r h r 2 ,r 3 )' 
^°(ri,r 2 ,r 3 ) 
\VZ s (r l ,r 2 ,r 3 ) / 



( cos 2 (f ) sin 2 (f) -sin(f) 
sin 2 (f) cos 2 (f) + sin(^n 



3 sln lf ) 



■ sin 



(¥) 



cos 



3 ' 
2k_ - 

3 ' 



■r 3 ,-r 2 ,-ri 



v$(-T 3 ,-T 2 ,- ri : 

V^(- r 3,-r2,-ri)y 



• (3-4) 



Similar relations hold, e.g. for the interaction vertices in the pNRQCD Lagrangian appearing at 
higher order in the multipole expansion in the relative vectors r,- (i =1,2, 3). 

The symmetry properties of the three heavy-quark system under exchange of the heavy-quark 
fields thus have a deep impact on the structure of the pNRQCD Lagrangian. They in particular in- 
duce relations between different matching coefficients in the effective theory and thereby constrain 
their form. 



4. The singlet static energy up to order ce In a s 



Besides studying the symmetries of the three heavy-quark system, we have used the effective 
field theory framework of pNRQCD to determine the leading ultrasoft contribution to the singlet 
static energy, which is of o^lna s , and to the singlet static potential, which is of order a^ln/x. Here, 



the symmetry relations (3.3) and (3.4) have served as an important check of the obtained result. 

Adding the newly determined leading ultrasoft corrections to the singlet static potential V s , 
known analytically at next-to-next-to-leading order (NNLO) [f§], the singlet static energy is now 
known up to order a s 4 In a s and reads 



£>i,r 2 ,r 3 ) 

«s 4 , 

In a s 

371 



V, 



NNLo( r b r 2) r 3 

(r 2 + r 3 ) 2 



l r i| 



+ 



l r 2| 



+ 



l r 3| 



1 |ri| + |r 2 | + |r 3 | 
4 |ri||r 2 ||r 3 | 
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1 

l r ll 



+ ■ 



l r 2| 



+ ■ 



(r2 + r 3 J 



l r i! 



+ 



l 



l r 2| 
1 

l r l| 



+ ■ 



l r 3| 



| 5 |ri| + |r 2 | + |r 3 | 
| 2 4 |ri||r 2 ||r 3 | 



1 



-ri • (r 2 + r 3 J 



1 



r i 



+ ■ 



1 



+ ■ 



1 



2|r 2 | 2|r 3 | 



; 5 |ri| + |r 2 | + |r 3 | 
4 |ri||r 2 ||r 3 | 



(4.1) 



l r 2| l r 3| 

1 1 

l r 2| |r 3 | 

In contrast to the static energy, the singlet static potential explicitly depends on the factor- 
ization scale /i separating soft from ultrasoft contributions. It is now known up to order af ln/i, 
where the quantity In /I is in general referred to as an ultrasoft logarithm. In a minimal subtraction 
scheme, it is given by 



V s (r u r 2 ,r 3 ;n) 



371 



In /I 



^NNLO ) 1*2 ) T 3 

(r 2 + r 3 ) 2 



+ 



l r il 



+ 



l r 2| 



+ ■ 



l r 3| 



1 |ri| + |r 2 | + |r 3 | 
4 



l r i|| r 2||r 3 | 



1 1 

— r + i— i 

r l r 2 



+ 



1 



(i*2 + r 3 ; 



r ll 



+ 



1 



l r 2| 

1 

r ll 



+ ■ 



l r 3 



r 3 



| 5 |ri| + |r 2 | + |r 3 | 
2 4 |ri||r 2 ||r 3 | 



1 



+ri -(r 2 + r 3 ) 



1 



+ ■ 



1 



r 2 



+ 



1 



r 3 



2|r 2 | 2|r 3 | 



5 |ri| + |r 2 | + |r 3 | 

2 + 4 



r i r 2 r 3 



(4-2) 



1 1 

|r 2 | |r 3 | 

Specializing to an equilateral geometiy, characterized by the single length scale r = |ri| = |r 2 | = 
|r 3 |, we have moreover managed to resum the leading ultrasoft logarithms that start appearing in 
the static potential at NNNLO to all orders by solving the corresponding renormalization group 
equations. For the singlet static potential this results in 



NNLO 



9 « s 3 (lA) ln «s(lA) 
/3 r a s (ju) 



(4.3) 



where j8o = 11 — 2/3«/, with n\ the number of light-quark flavors. Equation (4.3) provides the 
complete expression of the singlet static potential at next-to-next-to-leading-logarithmic (NNLL) 
accuracy in an equilateral geometiy. Corresponding results for the octet and decuplet potentials 
can be found in []2p. 
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